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ABSTRACT 

Recent high resolution numerical simulations have suggested that the in- 
terstellar atomic hydrogen clouds have a complex two-phase structure. Since 
molecular clouds form through the contraction of HI gas, the question arises as 
to whether this structure is maintained in the molecular phase or not. Here 
we investigate whether the warm neutral atomic hydrogen (WNM) can exist in 
molecular clouds. We calculate how far a piece of WNM which is not heated by 
the UV photons could penetrate into the cloud, and find that in the absence of 
any heating it is unlikely that large fraction of WNM survives inside high pressure 
molecular clouds. We then consider two possible heating mechanisms, namely 
dissipation of turbulent energy and dissipation of MHD waves propagating in the 
WNM inside the cloud. We find that the second one is sufficient to allow the 
existence of WNM inside a molecular cloud of size ~ 1 pc having pressure equal 
to ~ 10 X PisM- This result suggests the possibility that channels of magnetised 
WMN may provide efficient energy injection for sustaining internal turbulence 
which otherwise decays in a crossing time. 

Subject headings: Hydrodynamics - Instabilities - Interstellar medium: kinemat- 
ics and dynamics - structure - clouds 



Introduction 



There is strong observational evidence that the neutral ISM consists of two distinct 
temperature regimes, the WNM and CNM (e.g. Heiles and Troland 2003). This is a natural 
result of most static (e.g. Field et al. 1969, Wolfire et al. 1995), time- dependent (Gerola et 
al. 1974), and magnetohydrodynamical (e.g. Gazol et al. 2001; de Avillez & Breitschwerdt 
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2005; Inutsuka et al. 2005) models of the ISM because of the temperature-dependence of 
the interstellar cooling function. Recent hydrodynamical numerical simulations with a high 
numerical resolution approaching or resolving the Field length (Hennebelle & Perault 1999, 
2000; Koyama & Inutsuka 2000, 2002, 2004; Audit & Hennebelle 2005; Heitsch et al. 2005; 
Vazquez-Semadcni et al. 2005; Hennebelle & Passot 2006) have shown that the structures 
of the interstellar atomic hydrogen is likely to be highly complex. In particular in these 
simulations, it is found that i) the structures of cold neutral medium (CNM) are locally 
in near thermal pressure equilibrium with the surrounding WNM even if the flow is very 
dynamical on large scales, ii) the CNM structures have a velocity dispersion close to the 
sound speed of the WNM and in) the two phases appear to be highly interwoven. This 
picture of a dynamical two phase model appears to be compatible with the observations 
supporting the existence of phases as well as the observations which finds large fraction of 
thermally unstable gas (Heiles 2001) and large pressure fluctuations (Jenkins & Jura 1983). 

Since molecular clouds are supposed to form by contraction of atomic gas, the question 
of whether this complexity and the two-phase behaviour may persist even in molecular clouds 
arises immediately. We believe that the high turbulence which takes place in molecular clouds 
reinforces this idea ; since molecular clouds are surrounded by an HI halo, the turbulence 
should almost inevitably induces some mixing between HI and molecular gas. 

From an observational point of view, we believe that the presence of warm gas inside 
molecular gas is difficult to prove or disprove unambiguously because of the surrounding HI 
halo. However it is well established that molecular clouds arc clumpy and have a low filling 
factor (see next section). Both facts are compatible with the idea that WNM may exist in 
those clouds. Based on an analytic prediction of magnetic waves excited during a magnetic 
interaction of two magnetised clumps (Chfford & Elmegreen 1983) such two-phase model of 
molecular cloud has been considered by Falgarone & Puget (1986) (see also Elmegreen & 
Combes 1992). In their model, the giant molecular clouds are constituted of an ensemble of 
self-gravitating magnetised clumps which are surrounded by low pressure warm neutral gas 
which is heated either by ultraviolet photons or by the friction between ions and neutrals. 

Here we investigate the possibihty that WNM may exist even inside dense molecular 
clumps at pressure significantly higher than the mean ISM pressure. For this purpose, 
we consider various mechanisms and show the ranges of physical conditions by order-of- 
magnitude calculations. In the second section we describe the physical processes which are 
important for the thermal equilibrium of WNM and estimate how far a piece of WNM may 
penetrate into a molecular cloud before it cools down. In the third section, we consider the 
dissipation of MHD waves as a possible heating source of the WNM inside molecular clouds, 
and we calculate the corresponding thermal equilibrium curve. We also discuss the effect of 
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turbulence. The final section concludes the paper. 



2. Basics 

2.1. Notations and Assumptions 

We consider a molecular cloud of typical length L, mass M, volume V, pressure P, 
density p, and velocity dispersion a. Wc assume that the cloud follows the Larson's law 
(Larson 1981), i.e., M = 100 Mq x {L/lpcf and 



where a* = 0.4 km/s. Thus the turbulent energy of the cloud is about 1/2 x M a^. The 
mean value of the number density is about 

n ~ 570cm-2 x ( . (2) 

We assume that two phases in near pressure equilibrium fill the cloud as in the atomic 
gas clouds; one cold and dense phase {Tq = 10 K, pc/ni = uq > 10^ cm~^ where m is the 
mean mass of the particle ), and the other warm and diffuse phase (Twnm ^ 10^ K, uwnm — 
Tic X Tc/TwNM > 1 cm~^), and we define the phase density contrast Vp = nc/nwNM ^ 1000. 
If / is the volume filling factor of the cold component, then the volume occupied by the 
warm phase is Vwnm — — f) ^V- With these definitions, we have 

J. _ {n — nwNM) ^ ~ "^wnm) 
(nc - hwnm) nc 

As mentioned in the introduction, the simulations of turbulent atomic hydrogen clouds have 
shown that the phases are deeply interwoven. If this structure is preserved during the 
transformation process from atomic gas to molecular gas, channels of WNM gas permeating 
the molecular cloud should be naturally produced and would remain as a natural outcome 
of the multiphase structure. 

We will further assume that the whole cloud is permeated by a magnetic field B having 
roughly the same intensity within the two phases. In molecular clouds the observed average 
magnetic intensity is about 20//G for n — 1000 cm~^ and for higher densities, B is roughly 
proportional to ^/n (Crutcher 1999). As we discuss below, we will restrict our attention to 
clouds having a fiUing factor of about | . In this case, the density of the cold phase is about 
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2n and the WNM density is 2n/rp > 2 cm~^. Assuming that the mean magnetic intensity is 
nearly the same in the warm phase, this indicates that the magnetic field intensity in WNM 
embedded in two-phase molecular clouds is 



the same in the cold and the warm component is in good agreement with observations 
in the diffuse ISM (Troland & Heiles 1986). The magnetic intensity obtained by direct 
measurements of the magnetic field in the diffuse ISM indicates values of 5 ~ 6/iG in the 
WNM and in the CNM (Troland & Heiles 1986, Heiles & Troland 2006) where the mean 
density is respectively about 0.5 cm~^ and 50 cm~^. These measurements suggest values 
of B* smaller than 15 //G which may indicate lower filling factor or that our assumption 
of uniform magnetic field in the cloud does not hold exactly. In the following, we adopt 
B* ~ 10/iG that seems in better agreement with the magnetic intensities measured in the 
diffuse ISM. We also assume that the fluctuating component of the mag netic field, ^{6B^), 
is equal to the mean magnetic field, Bq = (B) which corresponds to equipartition between 
the magnetic energy of the mean and the fiuctuating part of the magnetic field. In the 
following we refer to the energy of the fluctuating component of the magnetic fleld within 
the WNM as the magnetic wave energy. Figure 1 shows a schematic picture that illustrates 
the model. 

Finally we also assume that the cloud is embedded in a diffuse gas of standard WNM 
having a density n^^^M — O-^ cm~^, a magnetic fleld of uniform component B^^t ~ 6/xG, 
and a fluctuating component y^JSB^ ~ Bext- The cloud therefore receives a flux of mag- 
netic energy from the external interstellar medium. This magnetic energy is produced at 
large scales through various mechanisms, such as clump magnetic interaction (Clifford & 
Elmegreen 1983, Falgarone & Puget 1986), supernovae explosions or galactic differential 
rotation (e.g. MacLow & Klessen 2004). 

According to equations (2) and (3), the flUing factor of the cold molecular component 
in a 10-pc cloud is very small (say, less than 0.1). It is therefore very likely that in such 
a cloud, the external radiation is not very different in most of the cloud volume, from its 
value outside the cloud. Therefore the conditions are likely to be similar to the standard 
ISM considered by Wolfire et al. (1995) so that WNM can exist inside such a large complex 
if the pressure in the rarefied area is comparable to the ISM pressure. On the other hand, 
for a 1 pc cloud, either / ~ 0.5 or ric is higher than 10^ cm~^ indicating a high thermal 
pressure. In both cases the UV background is not intense enough to heat the warm phase 
(see Section 2.2). Therefore in about one coohng time, WNM cools down into cold gas unless 
it is heated by another source of energy. In the following section we describe the relevant 
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where B* < 20/\/2 ^ lAfiG. Note that the assumption of magnetic intensity being roughly 
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physical phenomena (degree of ionization and recombination time, heating due to cosmic 
rays, UV and soft X-rays, cooUng rate and coohng time) and investigate how far a piece of 
WNM which is not heated could penetrate inside molecular clouds. 

2.2. Cosmic Rays, UV, and Soft X-Rays 

Cosmic rays, far ultraviolet (FUV) photons and soft X-rays are the main ionization 
sources of the WNM. They also constitute the main heating sources of the standard WNM. 

The value of the ionization rate due to cosmic rays is not well constrained and could be 
possibly very inhomogeneous. In molecular clouds, cosmic ray total (primary and secondary) 
ionization rate of C = 2 ~ 7x 10"^''' (e.g. Goldsmith 2001, Wolfire et al. 1995) are usually 
considered. However in diffuse clouds, values 40 times higher have been proposed by McCall 
et al. (2003) for explaining the high Hj^ abundance (see also, Le Petit et al. 2004). Recently 
Padoan & Scale (2004) proposed that confinement of cosmic rays by self-generated MHD 
waves may naturally accounts for these large variations. Since no measurement of cosmic ray 
ionization rate in the WNM is available in the litterature and since we arc proposing a model 
in which cold molecular gas is embedded in warm HI medium, the cosmic ray ionization rate 
relevant for our study is highly uncertain. In the following, we therefore consider a standard 
ionization rate of C = 3 x 10~^^ s~^ and discuss the consequences of higher values. The 
corresponding heating rate is 10~^^ erg s~^ (Goldsmith 2001). As we will see later this is too 
small to contribute significantly to the heating of WNM. 

As shown in Wolfire et al. (1995), the photoelectric heating from small grains and PAHs 
due to FUV photons is the most important source of heating of the standard WNM and is 
about 1.0 X 10~^^eG'o erg s~^, where e is the photoelectric efficiency and is equal to about 
<0.1 and Go is the incident FUV field normalized to Habing (1968) 's estimate of the local 
interstellar value. In the case of the standard ISM this heating appears to be sufficient to heat 
WNM of density up to 1 cm^"'. In the case of molecular clouds, the UV field is reduced by 
dust extinction by a factor ~ cxp(— 1.8 x A^) where is the extinction in visual magnitude. 
For a 1-pc cloud of average gas density 10^ cm^'^, this gives roughly ~ 1 reducing the 
incident FUV flux by a factor ~ 6. Note that since in our model the cloud is a multiphase and 
clumpy object, the FUV penetration should be enhanced by the diffusion processes (Boisse 
1990) and the incident FUV flux larger than this value. In any case, however, the heating 
due to standard FUV field is not sufficient to permit the existence of WNM at pressure 
larger than the standard ISM pressure and therefore cannot permit the existence of WNM 
inside most molecular clouds. Thus in this paper we assume, for simplicity, that the heating 
due to UV radiation field is negligible. Note that this approximation corresponds to the 
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most pessimistic assumption regarding the existence of WNM inside molecular clouds, and 
may not be appropriate for a low pressure isolated molecular cloud or for a cloud in the 
neighborhood of a strong UV source. 

As shown in Wolfire et al. (1995) the ionization and heating due to soft X-rays becomes 
neghgible for column densities larger than lO^'^ cm~^. Such a column density corresponds to 
a length of less than 0.03 pc for the cold component which we assume surrounds the warm 
phase that we are considering. 

2.3. Ionization and Cooling Rate 

The recombination rate of electron onto proton is about 2.6 x lO'^^T^'^ '^Ue s~^, where 
T4 = T/IOOOOX and Ue is the electron densities. This leads to an ionization fraction of 

X ^ njn^^u - 10- [j^^Y^) (^) (^) , (5) 

which is about 10 times smaller than in the standard WNM. The other possible contributions 
are ionization of heavy elements (mainly carbon) due to FUV and ionization of atomic 
hydrogen due to soft X-rays. However, the former is negligible since the carbon abundance 
is about 3 x 10^"^ and the latter is also negligible as explained in Section 2.2. Before the 
piece of WNM penetrates into the cloud, its ionization fraction is the standard value of the 
ISM, i.e., X ~ 0.1. Once it enters inside the cloud, it takes about one recombination time. 
Tree — IMyr to reach the equilibrium value within the molecular cloud, i.e x ~ 10~^. 

We have calculated the cooling rate by considering the standard cooling mechanisms 
of the atomic gas that are described in Wolfire et al. (1995), namely emission of Lyman a, 
metastable lines (CII, 01) and fine-structure lines (CII, 01) taking the values from HoUenbach 
& McKee (1989) and Wolfire et al. (1995) with heavy elements abundances corresponding 
to the solar neighborhood. Note that the ionization degree of carbon is very uncertain since 
as discussed in the previous section the UV flux is difficult to estimate. Therefore we make 
the assumption that carbon is fully ionized. This again corresponds to the most pessimistic 
assumption regarding existence of WNM inside molecular clouds since lower value of C"*" 
abundances will reduce the cooling (by about 20%). Note also that in spite of the fact 
that the whole cloud we are considering has a high opacity, the WNM embedded inside this 
cloud has a low column density and remains optically thin. Thus, the Lyman alpha photons 
emitted in the WNM propagate until they are absorbed by the dust into a piece of cold and 
dense gas. The corresponding energy is then reradiated by the cold component and finally 
leaves the cloud. Opacity effects are neghgible for the [01] 63 iim photons emitted by the 
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WNM and these photons escape the cloud. For the [CII] 158 //m photons, the opacity effects 
may depend on the exact structure of the cloud but remain modest (Tielens & HoUenbach 
1985). 

The main difference from the standard ISM calculations is about 10 times smaller ion- 
ization degree of the WNM. Therefore the radiative cooling is a few times smaller as well 
(see e.g. Dalgarno & McCray 1972), leading to a cooling rate of about 4 x 10~^® erg cm^ s~^ 
at T = 8000 K. 

Figure 2 shows the coohng time of the WNM when it enters into the cloud as a function 
of temperature. The density is obtained by assuming pressure equilibrium. The largest 
temperature corresponds to the case of a piece of WNM which has been suddenly compressed 
(and therefore heated) to the cloud pressure. Because of the cooling and the absence of 
heating, the temperature of the fluid element decreases monotonically with time, and the 
relevant cooling timescale is given by the largest value. 

The cooling time strongly depends on the ionization degree. However since it is always 
smaller than about 1 Myr and since the recombination time is about 1 Myr, the ionization 
degree is likely to be given by the full line (except maybe for P = Pism)- 

2.4. Penetration Length of WNM in Absence of Heating 

As discussed previously, heating due to FUV and soft X-rays is not efficient enough to 
provide a significant heating deep inside molecular clouds. We therefore calculate how deep 
the piece of WNM may penetrate inside the cloud without being heated. The corresponding 
length is given by the product of the WNM velocity by the cooling time. 

In numerical simulations of HI flows (Koyama & Inutsuka 2002, Audit & Hennebelle 
2005), it is found that the WNM velocity is higher by a factor of a few than the velocity 
dispersion of the CNM. We therefore expect that Mwnm — (Across x <^c, where across is a factor 
of a few. This is consistent with the idea that the cloud is continuously swept up by the 
surrounding HI gas which has an internal velocity dispersion equal to a fraction of the WNM 
sound speed, Cwnm — 10 km/s. 

Assuming that the WNM velocity inside the cloud is about 1 km/s (note that for 1-pc 
cloud, the Larson's law gives a velocity dispersion of about 0.4 km/s), we find with Figure 2 
that the piece of WNM cannot penetrate into the cloud deeper than 0.5 ~ 1 pc for P = Pjsm 
and 0.1 pc for P = 10 x Pism- Thus we must conclude that in the absence of heating WNM 
can exist inside low pressure molecular clouds (P ~ Pism) of size L ~ 1 pc but cannot 
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penetrate significantly inside high pressure molecular clouds of size larger than 0.1 pc. 



3. Heating Rate due to Mechanical Energy Dissipation 

Here we estimate the heating rate of the warm neutral phase due to mechanical energy 
dissipation in the cloud. First we consider the dissipation of magnetic waves that propagate 
into the warm phase inside the molecular cloud. Similar calculations have been performed 
by Scalo (1977) in the case of dense molecular gas and by Ferriere et al. (1988) in the 
case of standard WNM. Both found that the dissipation of magnetic waves can provide 
a substantial heating rate. We then estimate the amount of energy available from the 
dissipation of turbulent motions which are observed in the dense gas of molecular clouds. 



3.1. Dissipation of Magnetohydrodynamical Waves in WNM 

3.1.1. Wave Dissipation and Wave Energy 

The magnetic energy of the waves per unit volume is {6B'^)/87i = {{6 B"^) / Bq) Bq / Sir , 
where Bq is the mean magnetic field, and 6B is the fluctuating part of the magnetic field. In 
the following we will assume Bq ~ which corresponds to energy equipartition between 

the magnetic energies of the mean field and the fluctuating component. Assuming that the 
kinematic energy of the waves is comparable to the magnetic energy (as is the case for pure 
Alfven waves), the wave energy, 

_ {6B^Bl 

and the total amount of energy per particles of WNM available from the waves is therefore, 
^wnm/A^wnm X {5B'^)/An ~ {{dB^) / B^)B^/4:Ti/nwNM where A^wnm is the number of WNM 
particles within the cloud. 

The propagation of the Alfven waves in a weakly ionized gas has been studied by Kulsrud 
& Pearce (1969). They show that for a wavelength smaller than the critical wavelength, 

Acrit = ) (7) 

TdaPi 

the wave does not propagate, whereas for a wavelength larger than the critical wavelength, 
the wave propagates but dissipates in a timescale given by 

'"~^;i(27r/A)2' 
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where A is the wavelength, 7da is the friction coefficient between ions and neutrals, pn and 
Pi are the neutral and ion densities, and is the Alfven velocity. The expression of 
for a weakly ionized hydrogen gas has been derived by various authors (e.g. Osterbrock 
1961, Drainc 1980, Mouschovias & Paleologou 1981). Here we adopt the most recent value 
obtained by Glassgold et al. (2005), 7da = 5.7 x 10^"^ cm^ s~^ g~^(t;rms/kms~^)''''^^ where 
v-rms = ^/8kT/7rm is the mean thermal speed. This leads to 7da — 3.4 x 10^^ cm^ g~^ 
s-i(T/8000K)0-375 . 

Using the expression of the magnetic field (eq. [4]) and the ionization rate (eq.[5]), we 
obtain 

... . 0.0.,. (i^)- (^)-"- {^y^ (^) . „ 



3.1.2. Replenishment of Mechanical Energy 

Another important issue is on the mechanism for maintaining turbulent motions in 
molecular clouds. In the model that we are investigating, the WNM is heated by the dis- 
sipation of MHD waves. Therefore these waves must be continuously replenished from the 
outside, otherwise the heating source will disappear in the dissipation timescale. In the case 
of a single phase molecular cloud, Nakano (1998) pointed out that the dissipation timescale 
is shorter than the crossing time, which implies that internal turbulence cannot be easily 
driven from the outside. Here we examine three conditions that must be satisfied in order 
to enable energy replenishment from the outside. 

The first condition is that the timescale for the waves traveling into the WNM to reach 
the cloud centre, Across = {L/2)/va, must be smaller than or comparable to their dissipation 
timescale. The ratio between crossing timescale and dissipation timescale is given by 



Across _ Tt'^VaL _ TlXcritL 

Thus, the condition tcross/^da < 1 leads to A > Across where 
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Across = (7'"-^Acrit) ^ (11) 

~ 27 r^^^ f ^ ) 

~ ■ ^^llcm-^J VSOOOK; Vlpc/ \IQ p^G ) ' 

Note that this expression is valid only if L > Across- Therefore the smallest value of Across is 
obtained when L ~ Across and is about TrAcrit- 
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If the cloud is very anisotropic then the length of the shortest axis should be used in 
the expression of Across- 

The second condition is that the energy flux due to MHD waves which penetrate inside 
the cloud must be larger than the total energy dissipated by unit of time inside the cloud. 
This flux of energy cannot be easily computed but the maximum value can be estimated as 

Fe ^ 4:71(1 — /)^/^(L/2)^'t;A*-E'^ale where v*^^ and -E^avc ^^c respectively the Alfven velocity 
and the total energy of the waves outside the cloud and 47r(L/2)^ is the surface of the 
cloud. The factor (1 — Z)^^"^ takes into account the fact that only a fraction of the surface is 
constituted by WNM. The number of WNM particles is V{1 — /) nwNM, therefore one must 
have r„a,veK;(l ~ /) ^WNM < Fe which leads to 

p ^ plim _ ^ 6-Bcxt ((^-Bext) ^-|^2) 

wave- wave nwNM^ (47r)3/2^nwNM\AwJM^ ' 

where £^^ave — ^-^cxt/^^r and where the factor (1 — /)^/^ has been dropped for simplicity since 
it is on the order of unity unless the value of / is close to unity. The external gas density 
and magnetic intensity outside the cloud presumably vary significantly from place to place, 
in particular if the cloud forms dynamically. To estimate the largest heating which can be 
possibly due to input of external energy we adopt standard ISM conditions, i?ext — 6/iG, 
{SB^y.^) ~ B^^^^ and n^^M — 0.5 cm~^ which leads to 

r-, ^ 8 X 10-- ergs- ( ( f^V' ( -^p^V'' ( -^V , (IS) 



Bext J \6/iG/ Vlcm 3/ \0.5cm J \lpc^ 

For clouds which follows the Larson's law (eq. [2]), one finds that (^wnm/I cm-^)(L/lpc) ~ 
1, so that r|j,™yg ~ 8 X IQ-^^ ergs"^ is independent of gas density and cloud length. As we will 
see in following sections (see Fig. 4 and 5), this is sufficient to heat WNM inside molecular 
cloud of size L > 0.5 pc at the largest equilibrium pressures that are computed below. 

The third condition is related to the existence of the channels of magnetised warm gas 
which must permeate the molecular cloud in order for the energy to be replenished from the 
outside. Because of heat conduction the tunnels of WNM cannot be infinitely thin and in 
any case should be small compared to the size of the cloud. The smallest size of a piece of 
WNM embedded into cold gas is given by the Field length of WNM which is the typical size 
of a thermal front between the two phases. For smaller size the heat fiux due to thermal 
conductivity between the two phases will cool the WNM. The Field length is about 



~ V r ' ^^^^ 

where k{T) is the thermal conductivity, k{T) = 5/3C^rj{T) and 77 = 5.7 x 10-^{T/lKy/^ 
g cm"^ s"^ and F is the heating rate. In the standard WNM, the Field length is about 
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0.1 pc. As we will see in the next section, for the case that we are studying the heating 
rate can be more than 10 times higher than the heating rate of the WNM in standard ISM 
(see section 2.2) leading to a Field length smaller than 0.1/ vTo = 0.03 pc. The value is 
therefore small compared to the sizes of molecular clouds that enable the existence of the 
channels. Finally we also estimate the size of the smallest channels that can exist in spite 
of the thermal diffusivity during one dynamical time. If A^fp is the mean free path, then 
a particle undergoes a collision in a time of about Amfp/C's. During the collision with a 
colder particles, the WNM particle losses its thermal energy. Therefore the time required 
in order for the whole finger of WNM to cool down by collision with colder particles is 
(i?/Ainfp)^ X Amfp/Cg = -R^/AmfpCg whcrc R is the radius of the finger. If we require that the 
finger exists during a cloud crossing time, L/Cg, we find that R ~ y^LAmfp. For a density of 
n ~ 1 cm^^, Ajnfp — 10^^ cm. Therefore R ~ 0.03 pc for L = 1 pc which is thus comparable 
of the Field length within the cloud. 



3.1.3. Heating Rate 

Here we calculate the heating rate of the WNM due to MHD wave dissipation. Since 
tda depends on A, a wave powerspectrum has to be considered. 

We assume that the (isotropic) power spectrum of MHD waves E{k) has significant 
values only between /cmin = STr/Amax and /Cmax = STr/Amin, thus, the wave energy E'^ave — 
rfcmax £./^w^^ where we simply require 

Acrit ^ Amin < Amax ^ L. (15) 

The value of Amin is uncertain. On one hand, equation (11) gives the smallest wavelength 
which can be injected from the outside. On the other hand, wave steepening of long wave- 
length waves could generate waves of wavelengths close to Acrit- For the sake of simplicity 
we ignore this last possibility and assume that /cmax = Stt/ Across- As we discuss below, larger 
values of fcmax lead to larger heating. 

Since power spectra are often observed to be power-law, wc further assume S{k) = 
^o{kmm/ky where we may assume 1 < p < 2 according to the observations of molecular 
clouds. Note that the so-called Kolmogorov spectrum corresponds to p = 5/3. We have 



= r^^ Soik^n/kfdk = ^ [ 1 - ' ] . (16) 



-S'wave 



Obviously the energy is dominated by the waves of the largest wavelengths if p > 1. 
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The energy dissipation rate (per volume and time) is 

£^wave= / T\Adk= / So{k^n/ky,^k^dk. (17) 

Strictly speaking, depends on /c. However taking this into account would add more 
complexity to the final result without improving significantly its physical meaning. We 
therefore ignore this and write with equation (16) and using the expression for t^a 

P _ p /"''"""/'i. . /JAPh'^rJh - y -^wave^'i^min (femax/fcmin)^ ^ - 1 

^wave - ^0 2^^^^^ J^^^^ (fcm.n//c) k dk - ^ _ ^ X ^^^^^^ ^ _ ^f^^j ^^^y., ■ (18) 

Therefore in the case of the wave powerspectrum, the heating rate due to MHD wave dissi- 
pation can be written as 

'^max 

Twave — -S'wave — fi.Pi T )r'wave,o('^max) (1^) 

"'min 

where 

f{p, q) = ^ — X ^ — (20) 

rwave,o(A) is the heating that would be obtained with monochromatic waves of wavelength 
equal to A. Taking into account the expression of the magnetic field and the ionization rate, 
this leads to 



Vicm-3y Uoook/ 



rwave,o(A) ~ 8 X 10-2^erg s-^ (22) 

-1/2 f T / c y'^' / A / B* y nsB^y 
V3- 10-1^1 ) Vipc/ vio/^gJ \ W 

It is a significant heating which is significantly larger than the heating due to FUV in standard 
conditions. 

Figure 3 displays the numerical factor /(p, q) for typical values of p and q. Note that 
f{p,l) = 1 irrespective of the value of p. For a 1-pc cloud, the typical value of q is about 
L/Across — 5. This indicates that the actual value of F^ave is about 3 ~ 4 times F^ave.o- 
Note that integration upto /Cmax = 27r/Acrit instead of 27r/ Across; would lead to g ~ 40 and 
consequently to much larger heating that would violate the condition stated by equation 
(12). Since equation (19) appears to be complex, two asymptotic cases are discussed in the 
appendix. 
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Figure 4 shows the heating rate, F^ave in equation (19) for various density. For simphcity 
we set p = 5/3. Note that we have varied the slope of the power law, p between 1 and 3 
and found small departures from the heating obtained with the Kolmogorov spectrum as it 
is clear from Figure 3 for q < 10. The full lines arc for clouds following Larson's law so that 
n oc 1/L, while the dotted and dashed lines are for clouds of gas densities 1 and 3 cm~^, 
respectively. As expected the heating strongly depends on the cloud size L. 

3.1.4. Thermal Equilibrium 

We now calculate the thermal equilibrium for the WNM. We solve the energy equation 
nwNM^wavc = 'f^wNM^i'^) whcrc A(T) is the coohng function described in Sect. 2. 3 and 
Fwavc is given by equation (19) for various cloud sizes, namely 0.25, 0.5, 1 and 2 pc. Figure 
5 shows the thermal equilibrium curves. The points which have positive slope correspond to 
thermally stable states whereas those having a negative slope are unstable. The existence 
of WNM requires that the largest pressure reached by the thermal equihbrium curve should 
be larger than the pressure inside the molecular cloud. 

As expected the largest possible density at which WNM can exist depends strongly on 
the cloud size L. For L = 1 pc, pressures up to 10^ K cm^^ and densities up to ~10 cm~^ 
can be reached. For clouds of size L ~ 0.25 pc, which is close to the smallest wavelength 
TT X Acrit allowing energy replenishment, pressures up to 3x10^ K cm~^ and densities of about 
100 cm~^ can be obtained. We note however that, as shown in Figure 4, for L < 0.5 pc the 
heating is larger than T^^^^ which means that the corresponding energy is not available in 
standard ISM conditions. 

The exact value of the highest possible WNM pressure is more or less proportional to 
the heating which depends on few not-well-constrained parameters. In particular the amount 
of waves energy proportional to {SB'^) is a major source of uncertainties. However it should 
be noted that even with values ten times smaller, the highest WNM pressure would still be 
a few times the standard ISM pressure. 

3.2. Dissipation of Turbulent Energy 

Since molecular clouds arc turbulent it is unavoidable that turbulent energy is continu- 
ously dissipated. Here we estimate the heating rate of the WNM due to the turbulent energy 
dissipation. We note that the dissipation of turbulent energy in intermittent regions of energy 
dissipation within CNM and its consequences on the heating of the gas and the formation 
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of molecules has been investigated by Joulain et al. (1998). They find that temperatures 
up to 1000 K can be obtained in the dissipative structures. Padoan et al. (2000) explore 
the heating by ion-neutral friction in three-dimensional simulations of turbulent magnetised 
molecular clouds and find that this heating depends on the position and can be much higher 
than the heating due to cosmic rays. Recently warm molecular gas (rotational excitation 
temperature of 276 K) presumably heated by the dissipation of turbulent energy, has been 
observed by Falgarone et al. (2005). 

The mean rate of turbulent energy dissipation within the whole molecular cloud is 
estimated to be Ma'^/L. i.e., the available kinetic energy divided by the crossing time of the 
dense material. It is likely that most of this energy is dissipated and radiated away in the 
dense material, and only a fraction, r/turb; oi it can be used to heat the warm phase. We 
therefore estimate that the warm gas within the cloud can receive an energy rate due to the 
decay of the kinematic energy of rjturhMa^/ L. The ratio between the mass of the dense and 
warm components being given by ~ //(I — /) x Tp, the heating rate per particle in the WNM 
due to kinetic energy dissipation in molecular clouds is about //(I — /) x ?7turbrturb, where 
Tturb = rpTTia^/L. Using the values defined in Section 2.1, we obtain: 

r.„. = 10-ergs-.(^)'(A)"^ (23) 

This indicates that for a 1-pc cloud that follows the Larson's laws, the heating rate of WNM 
due to turbulence is small compared to the heating due to mhd waves dissipation. Tturh 
is displayed in Fig. 5 (dot-dashed thick line) which shows that the turbulent heating may 
become dominant for L > 5 pc. However rjtuvb is likely to be small since dissipation occurs 
within the cold component which carries the turbulent energy. We therefore expect that 
Vturh^turb IS significantly smaller than r^,,avc- It is therefore unclear whether the turbulent 
energy dissipation is important to heat the WNM. It is however possible that, as proposed 
by Clifford & Elmegreen (1983) and Falgarone & Puget (1986), in a giant molecular cloud 
the kinetic energy of the translational motions of cold clumps could be the dominant energy. 
In this case, the motion of the cold clumps excites Alfven waves which then cascade and 
dissipate providing heating of the WNM at smaller scales. 

Another difficulty with the turbulent energy is that it dissipates within one cloud cross- 
ing time (e.g., Mac Low & Klcsscn 2004) and needs to be continuously maintained by some 
forcing which is not identified yet. If the WNM exists deep inside molecular clouds as we 
propose here, the high Alfven velocity in the warm phase enables the injection of the MHD- 
wave energy from the outside of the molecular clouds, and it might be transmitted to the 
cold component of molecular clouds contributing to the turbulent motions. Since both en- 
ergy densities have comparable values, E'wave — -E'turb) ^-nd since the fiux of magnetic energy 
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is larger than the flux of turbulent energy, E'wave x K,wnm ^ -Eturb x cr, the energy could 
in principle be injected inside the cloud through MHD waves propagating inside WNM and 
then partly transmitted through MHD interaction to the cold component leading to near 
equipartition of the energies. 

4. Discussion and Conclusion 

We have investigated the possibility that WNM may exist in molecular clouds by consid- 
ering the cooling time of the WNM fluid particle and the heating rate due to the mechanical 
energy dissipation into the cloud. 

Our estimate indicates that the WNM when it enters into a high pressure molecular 
cloud (P ~ 10 X Pism) cools too rapidly to allow the existence of the warm phase inside a 
cloud of size larger than ~ 0.1 pc unless it is heated by some process. On the contrary, in 
a low pressure cloud (P ^ Pism) the WNM can penetrate deep into a cloud of size up to ~ 
1 pc. The dissipation of the turbulent energy of the cold component of the molecular cloud 
and the photoelectric heating from small grains and PAH may provide enough heating only 
for low pressure molecular clouds and unlikely for molecular clouds having pressure 10 times 
the pressure of the ISM unless the cloud is very turbulent (velocity dispersion about 2 times 
the value of the Larson's law) or located near a strong UV sources {Gq ~ 10). 

In contrast, the dissipation of the MHD waves seems to be a promising mechanism to 
maintain warm atomic gas inside high pressure molecular clouds. The value of the highest 
pressure at which WNM may exist depends on the cloud size, density, magnetic intensity 
and ionization degree. We flnd that for a 1-pc cloud, WNM may exist up to pressures of 
about 10 X Pism- 

The main difference with the previous estimate for the case of the standard ISM (e.g. 
Ferriere et al. 1988) is that, since the ionization is about 10 times smaller than in the standard 
ISM, the coohng rate is few times smaller, whereas the heating rate is larger because the 
dissipation time is smaller. Also the magnetic fleld is stronger and the wavelengths smaller. 

Finally we suggest that if channels of warm diffuse and magnetised gas do exist and 
permeate the molecular clouds, then the dissipated energy can be replenished from the 
outside because of the high Alfven velocity in the WNM. We speculate that such channels 
may also help to inject energy continuously into the dense component of the molecular cloud 
helping to sustain the turbulent motions that otherwise decay in a crossing time (e.g., Mac 
Low & Klessen 2004). 
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A. Asymptotic Limits for Heating due to MHD Wave Dissipation 

Since equation (19) appears to be complex two asymptotic cases are discussed. 

The first asymptotic case is when the size of the cloud, L, is close to Amin = 27r x Acrit 
then /cmax — 27r/Amin and /cmax = tt^min whcrc a is a numerical factor of a few, and we have 

1 a4/3 _ 1 
where r™^^ = rwave,o(A = Amin) and is given by 

.4 X 10-3,,^,-. (^)- (_r_y- (_^)"^ (_£_)\ 

Since r„avc,o decreases with A, r™!^^ is an upper limit for the heating rate due to MHD 
dissipation. Unlike F^g^^g, F™^g increases with the cosmic rays ionization rate and with the 
gas density. This is due to the fact that Amin decreases with these two parameters and that 
the dissipation time is oc A~^. For a = 2, we have Fwave — 0.5 x F^l^^. 

The second asymptotic case that we consider is when the size of the cloud is much larger 
than Amin, ^max = 27r/Across ^ ^min = 27r/L. Then we can write 

T-wave — rwave,o(A = Across) ^ 2 ( ^ j ' i-^^) 

Gombining this equation with equation 11, we obtain 

J, \ -0.242 / ^ \ '^/^ f L \ f B* \ ^^^^ 



r ~ 1 6-1 n 

'- wave — -'-■'J J-U 



(A4) 
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This expression depends only weakly on the gas density and ionization rate. It decreases with 
the size of the cloud, L, less rapidly than in the monochromatic case stated by equation (23). 
Note that for /cmin = 27r/Ainin instead of k^in = 27r/Across; one obtains F^ave oc 't-wnmC^^^-^c 
indicating an even shallower dependence on the cloud size. 
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Fig. 1. — Schematic picture illustrating the model of multiphase magnetised molecular 
clouds. 
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Fig. 2. — Cooling time of the WNM when it enters a molecular cloud of pressure equal to 
the ISM pressure, Pism (3 top curves), and 10 x Pism (3 bottom curves). Full line is for an 
ionization degree, a;, of 0.1, dotted line is for x =0.01 and dashed line is when the ionization 
equilibrium is assumed. 
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3. — The numerical factor f{p, q) in equation (19). 
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Fig. 4. — Heating rates due to MHD waves dissipation F^ave stated by equation (19). Full 
lines are for clouds following Larson's laws so that n cc 1/L. Dotted and dashed lines are for 
clouds of gas densities 1 and 3 cm^"^ respectively. The dot-dashed line displays the heating 
due to the dissipation of turbulence. The triple dot-dashed line displays the heating that 
would be obtained by dissipating all the available magnetic wave energy in standard ISM 
conditions (see eq. [12]) for clouds following the Larson's laws. 
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Fig. 5. — Thermal equilibrium curves for a heating rate due to MHD wave dissipation r^ave 
and for four cloud sizes. Only densities corresponding to the warm phase are shown. The 
points which have positive slope correspond to thermally stable state whereas the one having 
negative slope are unstable. The existence of WNM requires that the pressure is lower than 
the largest pressure reached by the thermal equilibrium curve. 



